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ABSTRACT 

We have calculated fast direct spectral model fits to two early-time spectra of the 
Type-II plateau SN 1999em, using the SYNOW synthetic spectrum code. The first is 
an extremely early blue optical spectrum and the second a combined HST and optical 
spectrum obtained one week later. Spectroscopically this supernova appears to be 
a normal Type II and these fits are in excellent agreement with the observed spectra. 
Our direct analysis suggests the presence of enhanced nitrogen. We have further studied 
these spectra with the full NLTE general model atmosphere code PHOENIX. While we 
do not find confirmation for enhanced nitrogen (nor do we rule it out), we do require 
enhanced helium. An even more intriguing possible line identification is complicated 
Balmer and He I lines, which we show falls naturally out of the detailed calculations 
with a shallow density gradient. We also show that very early spectra such as those 
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presented here combined with sophisticated spectral modeling allows an independent 
estimate of the total reddening to the supernova, since when the spectrum is very blue, 
dereddening leads to changes in the blue flux that cannot be reproduced by altering the 
"temperature" of the emitted radiation. These results are extremely encouraging since 
they imply that detailed modeling of early spectra can shed light on both the abundances 
and total extinction of SNe II, the latter improving their utility and reliability as distance 
indicators. 

Subject headings: radiative transfer — stars: atmospheres — supernovae: SN 1999em 

1. Introduction 

SN 1999em was discovered by the Lick Observatory Supernova Search (IAUC 7294) on Oct 
29.44 UT and confirmed by the BAO supernova group on Oct 29.7 UT. It was followed spectro- 
scopically and photometrically by ground-based observatories and was the subject of a Supernova 
Intensive Studies (SINS) Hubble Space Telescope (HST) observation on Nov 5, 1999. SN 1999em 
is in the SABc galaxy NGC 1637, with a heliocentric recession velocity of 717 km s _1 . The first 
spectrum taken on Oct 29, 1999 is extremely blue and displays strong Balmer features typical of 
Type II supernovae with a characteristic velocity of around 10,000 km s _1 . 

SN 1999em appears to have been a normal Type IIP, likely the result of a core-collapse in a 
star with a massive hydrogen envelope (Leonard et al. 2000). The early and weak radio emission 
(IAU Circulars 7318 and 7336) hints that the star did not experience significant mass loss shortly 
before the explosion, which is also thought to be the normal situation for SNe IIP (see Weiler et al. 
2000, and references therein). The lack of circumstellar interaction makes SNe IIP "clean" cases 
to model, resulting in more reliable atmospheric modeling which in turn makes SNe IIP ideal for 
cosmological use. Because SN 1999em was observed earlier, and over a wider wavelength range 
than any other normal bright SN IIP, it is a perfect test bed for our understanding of SNe IIP. 
These supernovae are of particular interest as distance indicators. 

2. Observations 

The observed spectra we use in this paper were obtained on Nov 4, 1999 at the Fred L. Whipple 
Observatory (FLWO), and on Oct. 29, 1999 at CTIO. Leonard et al. (2000) find that on Oct. 29 
the V light curve is rising. The point at Nov. 4-5 is on a broad maximium prior to the plateau. 
The light curve settles down to a plateau phase on about Nov. 20. 
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3. SYNOW Models 

We have used the fast, parameterized supernova synthetic spectra code SYNOW to make an 
initial investigation of line identifications and expansion velocities. The code is discussed in detail 
by Fisher (2000) and recent applications include Fisher et al. (1999), Millard et al. (1999), and 
Hatano et al. (1999a). SYNOW spectra consist of resonant-scattering profiles superimposed on 
a blackbody continuum. For the two synthetic spectra presented here, the radial dependence of 
all line optical depths is a power law of index 8. The observed spectra have been corrected for a 
redshift of 717 km s _1 and for reddening of E(B — V) = 0.05 mag (see below). 

In Figure 1 the CTIO observed spectrum of October 29 is compared with a synthetic spec- 
trum for which the continuum blackbody temperature and the excitation temperature T exc are 
13,000 K, the velocity at the photosphere V p h ot is 11,000 km s^ 1 , and only lines of hydrogen, He I, 
and N II are considered. The relative strengths of the lines of an ion are fixed by local thermody- 
namic equilibrium (LTE) at the excitation temperature, but the absolute strengths of the lines of 
each ion are controlled by a free parameter. The strength of the hydrogen lines has been chosen to 
make H/3 and H7 about the same as in the observed spectrum; the relative optical depths of the 
hydrogen lines are of course fixed by atomic physics so the poor fit to Ha is due to a deficiency of 
the resonant-scattering source function for this transition. The only noticeable line of He I in the 
synthetic spectrum is that of A5876. Lines of N II have been introduced so that A4623 and A5679 
can account for the absorptions near 4500 and 5500 A, in which case N II A5029 also affects the 
synthetic spectrum slightly near 5000 A. (In SYNOW when a particular species in introduced, all 
the lines of that species are included with their relative strengths set by assuming LTE populations, 
but we focus here on the strong optical features which are relevant to the observed spectra.) The 
hydrogen and He I identifications are definite, and (in anything remotely like LTE) we can offer no 
alternative to the N II lines. However, see § 4.1.1 regarding the possibility that some of the "N II" 
lines may actually be due to hydrogen. 

In Figure 2 the HST observed spectrum of November 5 is compared with a synthetic spectrum 
that has T bb = T exc = 9500 K, v phot = 8000 km s _1 , and lines of hydrogen, Ca II, Mg II, Fe II, and 
Ni II. The resonance-scattering approximation gives good fits to the features produced by H/3, H7, 
and H<5. Ca II contributes only the blend due to A3945 (the H and K lines) and Mg II contributes 
little other than the blend due to A2798 (the h and k lines). All other features in the synthetic 
spectrum are produced by Fe II and, to a lesser extent, by Ni II. The hydrogen, Ca II, and Mg II 
identifications are definite. The ultraviolet spectrum evidently is mainly a blend of singly-ionized 
iron-peak lines; Fe II definitely is present and Ni II probably contributes significantly. Other 
iron-peak ions that have not been introduced here may also affect the observed spectrum. 
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4. PHOENIX Models 
4.1. Oct 29 

The Galactic extinction to NGC 1637 corresponds to a color excess of ^(.B — V) = 0.03 mag, 
but the observed spectrum of SN 1999em shows a distinct NaD interstellar absorption line that 
may indicate additional reddening in the parent galaxy. We have calculated a grid of detailed fully 
line-blanketed PHOENIX models in order to determine the "temperature" of the observed spectrum. 
With most supernova spectra we can trade off higher temperatures for larger reddening within 
prescribed limits (such that the atmosphere does not become so hot or cool that strong unobserved 
lines would predominate). Here, with this very early, very blue spectrum, we found that we were 
not able to freely exchange temperature for reddening, so that the reddening was well determined 
to be E(B — V) ~ 0.05 and the "temperature" was about 11,000 K. Figure 3 displays a fit with 
^rnodel = H>000 K and solar compositions; the observed spectrum has been dereddened with 
E(B — V) = 0.05 mag. The abundances were taken to be solar throughout the model atmosphere 
and the density was assumed to follow a power-law with p oc r~ n and n = 7, somewhat shallower 
than in the SYNOW fits, though we don't regard the difference as being too significant, since the 
true density structure is unlikely to follow an exact power-law. We also present models with n = 9 
(see Figure 8). 

The fit is rather good, and the ^ mo( jel cann °t be significantly reduced since the Ca II H+K 
line becomes very strong in the synthetic spectrum with ^ mo del = K, but it is weak in 

the observed spectrum. Also, the observed feature due to He I A5876 is not reproduced in the 
synthetic spectrum. Hotter models are not ruled out but at ^ moc | e i = 12, 000 K the Ca II H+K 
feature becomes weaker and there is a clear flux deficiency in the red. Normally one expects that 
one can trade off reddening for temperature in the red part of the spectrum, and this is indeed the 
case within a certain range of reddening. We assume a standard (Cardelli et al. 1989) reddening 
law with R = 3.1. However, since the spectrum depends strongly on the strength of the Ca II 
H+K feature, we find that ^ moc i e i = H, 000 ± 500 K; the T mo ^ e Y = 12, 000 K model requires an 
extinction E(B — V) ps 0.15 to reproduce the red flux, but it then does a poor job in the blue. 
Thus, we find that E(B - V) » 0.05 and E(B - V) < 0.15. We show a model with E(B - V) = 0.10 
below. 

We also found that we could not reproduce the He I A5876 feature by enhancing the amount of 
gamma-ray deposition. Therefore, guided by the indications of enhanced N in the SYNOW fits, we 
calculated a series of spectra where we enhanced He by a factor of 2.5 (reducing H appropriately), 
and we enhanced N by a factor of 10 (reducing C+O). This fit is shown in Figure 4, where now the 
He I A5876 fits well, but ^ mo( j e l = Hj 500 K in order to keep the Ca H+K feature from becoming 
too strong in the synthetic spectrum. Figure 5 with 2^ mo( j e i = 12, 000 K, the same compositions 
and E(B — V) = 0.10 shows that the effects of the Ca H+K feature can be reduced by increasing 
the model temperature and increasing the reddening modestly. Therefore we conclude that there 
is strong evidence for enhanced helium in SN 1999em [a similar result was found for SN 1987A 
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(Arnett et al. 1989; Eastman & Kirshner 1989; Lundqvist & Fransson 1996; Blinnikov et al. 2000)]. 
We note at these early times that the material is hot enough for the optical helium lines to be 
excited even in LTE (see Hatano et al. 1999b), so while gamma-ray deposition was needed in e.g. 
SN 1993J (see Baron et al. 1995a, and references therein), it is not required at these early times 
(nor is it precluded at later times, although we did not find it necessary for the spectra a week 
later) . 

4-1.1. Complicated Balmer and He I Lines 

We attempted to confirm the evidence in the SYNOW fits that nitrogen is also enhanced and 
that N II contributes to the observed features at A5679, A5029, andA4623. To do this we reduced 
the abundance of nitrogen from its solar value by a factor of ten. This synthetic spectrum showed 
that the A5679 and A5029 features do decrease, but the feature that is attributed to A4623 in the 
SYNOW fits is still present (this is not surprising, since the feature also appears in Figure 3 with 
solar compositions). Guided by the LTE line spectra presented by Hatano et al. (1999b), we reduced 
carbon by a factor of ten from its solar value and the 4600 A feature remained. Finally, reducing 
both carbon and oxygen by a factor of ten from their solar values (leaving N fixed) still produces 
a feature with if anything a better shape than in previous calculations. Thus, with the "usual 
suspects" the source of the 4600 A feature remains unidentified, and the evidence for enhanced 
nitrogen is weakened, but it does a good job on the A5679, A5029, features. 

The lack of a strong candidate for the 4600 A feature leads us to consider more exotic line 
identifications. Careful examination of Figure 3 reveals that even with solar compositions features 
appear in the synthetic spectrum at the approximately correct wavelength positions. This led us to 
test the hypothesis that the features are actually produced by "complicated P-Cygni" profiles, that 
is the usual wide P-Cygni profile with the peak centered near zero velocity and a second P-Cygni 
profile with a second absorption minimum around 20,000 km s _1 . Figure 6 displays the synthetic 
calculation with the metalicity reduced to Z = Zq/100. One can be confident that all the lines in 
this synthetic spectrum are produced solely by hydrogen and helium and it is clear that each line 
is associated with high-velocity absorption, which matches the observed lines very well. (It may be 
a true P-Cygni feature, with both a peak and a dip, but it is probably premature to decide that at 
this point.) This is due to the shallow density gradient and high model temperature which keeps 
the actual electron temperature above 8000 K all the way to the highest velocity in the model, 
37,500 km s _1 . Complicated NLTE effects vary the Balmer level populations in the mid-velocity 
range, producing the double P-Cygni lines, i.e. there exist two line forming regions for hydrogen 
and helium. This is the first time that we have encountered this result and it shows that particularly 
in differentially expanding flows, one cannot be certain of line identifications until careful synthetic 
spectral models have been calculated. Thus, we therefore believe that the features that we have 
heretofore attributed to N II (and which would require enhanced N) are in actuality due to H I and 
He I, which form secondary features at high velocity. Since this appears in our modeling without 
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any special parameters (i.e. Figure 3), we believe this is a true NLTE effect produced simply by 
the shallow density gradient. 

4.2. Nov 4 5 

Guided by our results in § 4.1 we have computed full NLTE models for the epoch observed with 
HST using both solar and enhanced He+N compositions (since we have not ruled out enhanced N, 
we will continue to use it for expediency). Figure 7 shows our best fit for solar compositions and 
n = 7 with r moc j e j = 8100 K. Figure 8 displays our best fit for the enhanced He+N compositions, 
where we have also found the fit is improved by setting T moc j e j = 9000 and n = 9. Also displayed 
is a ground-based optical spectrum obtained on Nov. 4, 1999 at the FLWO, which lends support 
for enhanced He. All the models show unobserved strong features in the range 2000-3500 A, which 
could be due to iron-peak elements. Thus, it is possible that the metallicity of SN 1999em was 
somewhat less than solar, but we will study this question in future work. 

5. Conclusions 

We have shown that both direct synthetic spectral fits and detailed NLTE models do a good job 
of reproducing the observed optical+UV spectra of SN 1999em. We have shown that detailed NLTE 
spectral modeling of very early spectra of Type II supernova can provide an independent estimate 
to the total reddening as well as abundance information. Specifically we find that E(B — V) ~ 
0.05 — 0.10 and E{B — V) <^ 0.15 mag, and we find strong evidence for helium enhanced by at least 
a factor of 2 over the solar value. We obtain the striking result that observed features are likely 
due to "complicated P-Cygni" lines of hydrogen and helium, and thus one needs to be especially 
careful with line identifications in hot differentially expanding flows. That the observed complicated 
profiles fall right out of the PHOENIX calculations is remarkable and lends strong support to the 
reliablity of the modeling. 

The fact that planned supernova search programs will preferentially find supernovae very 
early, combined with our results that the extinction can be determined by modeling the very early 
spectra, lends promise to the use of Type II supernovae as distance indicators, through the use of 
phenomenological applications like the "Expanding Photosphere Method" (Baade 1926; Kirshner & 
Kwan 1974; Branch et al. 1981; Schmidt et al. 1992; Eastman et al. 1996) or through sophisticated 
modeling of the "Spectral-fitting Expanding Atmosphere Method" (Baron et al. 1993, 1994, 1995a,b, 
1996). 
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Fig. 1. — SYNOW fit to Optical Oct. 29 spectrum. Observed spectra have been de-redshifted by 
717 km s _1 in this and subsequent figures. The observed spectrum has been de-reddened using 
E(B — V) = 0.05. 
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Fig. 2.— SYNOW fit to the UV+optical HST spectrum obtained on Nov. 5, with the optical 
spectrum obtained at the FLWO on Nov. 4. The HST and optical spectra overlap at around 
3600 A. The observed spectra has been de-reddened using E(B — V) = 0.05. The fit quality should 
be judged primarily in the optical. 
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Fig. 3. — PHOENIX fit (thick line) to the optical spectrum obtained on Oct 29 at CTIO. The model 
temperature ^ mo( j e l = 11,000 K and the compositions are solar. The red flux (A ^ 6500 A), 
indicates that the model temperature is about right. In this and subsequent PHOENIX plots the 
following species were treated in full NLTE: H I, He I-II, CNO I-III, Mg II, Ca II, Na I-III, Si I-III, 
S I-III, Fe I-III, Ni I-III, and Co I-III. 
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Fig. 4.— PHOENIX fit (thick line) to the optical spectrum obtained on Oct 29 at CTIO. The model 
temperature T mQ( ^ e i = 11500 K and the compositions are enhanced in helium and nitrogen. With 
this choice of composition and temperature the He I A5876 feature is well fit, and the Ca H+K 
feature is too strong, indicating that the model temperature is too low. The red flux indicates that 
the model temperature is too low; thus, a slightly higher reddening may reconcile the red and blue 
ends of the spectrum. The observed spectra have de-reddened using E(B — V) = 0.05. 
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Fig. 5.— PHOENIX fit (thick line) to the optical spectrum obtained on Oct 29 at CTIO. The 
model temperature is ^ moc iel = 12, 000 K and the compositions are the enhanced ones in Figure 4. 
The higher temperature of this model weakens the Ca II H+K feature seen in Figure 4, and the 
extinction has been taken to be E(B — V)= 0.10 which gives a reasonable fit to the red flux. 
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Fig. 6. — PHOENIX model (thick line) with Z = Z Q /100 compared to the optical spectrum obtained 
on Oct 29 at CTIO. All of the lines in the synthetic spectrum are produced solely by hydrogen and 
helium. The number abundances of the next most abundant species (C+O) are 10~ 6 smaller than 
those of hydrogen and helium. The reddening has been taken to be E(B — V) = 0.10. The inset 
shows a close up of H/3 and it is clear that there is a shoulder in both the observed and synthetic 
spectrum. 
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Fig. 7. — PHOENIX fit (dashed line) to the UV+optical HST spectrum obtained on Nov. 5, with 
the optical spectrum obtained at the FLWO on Nov. 4 attached (redward of 3600 A). The model 
temperature 2^ mo( j e l = 8100 K and the compositions are solar. The observed spectra have been 
de-reddened by E(B -V) = 0.05. 
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Fig. 8. — PHOENIX fit (dashed line) to the UV+optical HST spectrum obtained on Nov. 5. The 
model temperature 2^ m odel = 9000 K, n = 9, and the compositions have He and N enhanced. Also 
displayed is the optical spectrum (redward of 3600 A) obtained at the FLW Observatory on Nov 
4, which strengthens the case for enhanced helium. The observed spectra have been de-reddened 
using E(B — V) = 0.10. 



